Impossibility of superluminal travel in Lorentz violating theories 
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Warp drives are space-times allowing for superluminal travel. However, they are quantum me- 
chanically unstable because they produce a Hawking-like radiation which is blue shifted at their 
front wall without any bound. We reexamine this instability when local Lorentz invariance is vio- 
lated at ultrahigh energy by dispersion, as in some theories of quantum gravity. Interestingly, even 
though the ultraviolet divergence is now regulated, warp drives are still unstable. Moreover the type 
of instability is different whether one uses a subluminal or a superluminal dispersion relation. In the 
first case, a black-hole laser yields an exponential amplification of the emitted flux whereas, in the 
second, infrared effects produce a linear growth of that flux. These results suggest that chronology 
could still be protected when violating Lorentz invariance. 

PACS numbers: 04.62. +v, 04.60.Bc, 04.70.Dy 



I. INTRODUCTION 

Warp drives [1] allow, at least theoretically, to travel at 
arbitrary high superluminal speeds and consequently to 
travel in time [2]. However, besides the fact that they re- 
quire matter distributions violating positive energy con- 
ditions [3-6], they are quantum mechanically unstable 
because they possess a white hole horizon and a Cauchy 
horizon on which the renormalized stress-energy tensor 
blows up exponentially [7]. In this paper, we reexam- 
ine the question of their stability when postulating that 
Lorenz invariance is broken at ultrahigh energy. In this 
we were motivated by the fact that nonlinear dispersion 
relations remove Cauchy horizons and regulate the fluxes 
emitted by white holes [8]. 

A superluminal warp drive metric describes a bubble 
containing an almost flat region, moving at some constant 
speed vq > c within an asymptotically flat space-time: 



ds 2 = -c 2 dt 2 + [dx - v(r)dt] + dy 2 + dz 2 



(1) 



where r = y (x — vot) 2 + y 2 + z 2 is the distance from 
the center of the bubble. Here v = vof(r), with / a 
smooth function satisfying /(0) = 1 and f(r) — > for 
r — > oo. Along the direction of motion, the backward 
and forward loeii where v(r) — c behave, respectively, 
as a future (black) and past (white) event horizon [9]. 
In fact, the Hawking flux emitted by the black horizon 
accumulates on the white horizon while being unbound- 
lessly blueshifted. However, since the whole analysis rests 
on rclativistic quantum field theory, one should examine 
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whether the warp drive's instability is peculiar to the 
local Lorentz symmetry. Although current observations 
constrain to ultrahigh energy a possible breaking of that 
symmetry [10], one cannot exclude this possibility which 
has been suggested by theoretical investigations [11-13]. 

For the sake of simplicity we work in 1 + 1 dimensions 
and ignore the transverse directions y and z. Defining a 
new spatial coordinate X = x — vot, the metric becomes 



ds 2 = -c 2 dt 2 + [dX - V(X)dt] 



(2) 



where V(X) = vo(f(X) — 1) is negative. In this space- 
time, dt is a globally defined Killing vector field whose 
norm is given by c 2 — V 2 : it is timelike within the bubble, 
its norm vanishes on the two horizons, and it is spacelike 
outside. One thus gets three regions: L, C, and R (see 
Fig. 1), separated by two horizons xbh < < xwh, 
which are, respectively, the black and the white one. 

We now consider a massless scalar field with a quartic 
dispersion relation. In covariant terms, its action reads 



S4 



d 2 Xy 



A 2 



(3) 



where W = g^ v +u^u v is the spatial metric in the direc- 
tion orthogonal to a unit timelike vector field . This 
extra background field specifies the preferred frame used 
to implement the dispersion relation [14]. In the present 
settings should be given from the outset, while in con- 
densed matter the preferred frame is fixed by the fluid 
flow [15]. Inspired by this analogy, we choose u M to be 
(ljV") in the t,X frame, i.e. stationarity is preserved. 
Then the aether flow is geodesic and it is asymptotically 
at rest in the t, x frame of Eq. (1). The sign ± in Eq. (3) 
holds for superluminal and subluminal dispersion, respec- 
tively. Using Eq. (2) and u M = (1, V), the wave equation 
is 



(dt + d x v) (d t + vd x ) -d%± j^d x 



= 0. (4) 



*BH 



FIG. 1. Velocity profile for a right-going warp drive in the 
Painleve-Gullstrand [28] coordinates of Eq. (2). Two super- 
luminal asymptotic regions L and R are separated by a black 
and a white horizon from a compact internal subluminal re- 
gion C. The Killing field dt is spacelike in L and R, lightlike 
on both horizons, and timelike in C . 



Because of stationarity, the field can be decomposed in 
stationary modes <fi = ^ dui e~ lult 4> u , where oj is the con- 
served (Killing) frequency. Correspondingly, at fixed u) 
the dispersion relation reads 



(w - Vk u f = k; 



+ *Sl = o 2 
A 2 ~ ±: 



(5) 



where k u (X) is the spatial wave vector, and f2 the co- 
moving frequency, i.e. the frequency in the aether frame. 
The graphical solution of Eq. (5) is plotted in Fig. 2. For 
super luminal dispersion and |V| < 1, there are two real 
roots k%) describing left- and right-going waves 
</>"), and two complex ones (fc£, k^) describing a spa- 
tially growing and decaying mode 0^). For |V| > 1, 
there is a cutoff frequency w max [8] below which the com- 
plex roots turn into real ones (k^\ fc£ 2 ^) with negative 
VL. Correspondingly there exist two additional propagat- 
ing waves {4>^-J)*t (4>-l,)* w ^ n ne g a ti ve norm. When 
the dispersion relation is subluminal, the two extra roots 
correspond to positive norm modes that are trapped in 
the region with \V\ < 1. 




FIG. 2. (color online). Graphical solution of Eq. (5) for super 
(left panel), and subluminal dispersion (right panel). In both 
panels, the straight lines represent ui — Vk for \V\ < 1 (solid) 
and \V\ > 1 (dashed). The curved lines represents ±f2±(fe). 
On the left, red (blue) dots in the upper (lower) half plane 
refer to roots with positive (negative) 0+ which correspond 
to positive (negative) norm modes. 



FIG. 3. Asymptotic decomposition in plane waves ip. 
of the in mode 



(i),L/R 



Note that only <p_lj nas group 



(i) 



velocity directed toward the horizons, with wavevector k 



II. SUPERLUMINAL DISPERSION RELATION 

In a geometry with two infinite asymptotic "super- 
luminal" regions, for each u < w max , 4 asymptotically 
bounded modes [8] can be defined. Moreover, by exam- 
ining their asymptotic behavior, an in and an out bases 
can be defined by the standard procedure: each in mode 



i(0. 



(out mode 



) possesses a single asymptotic 



branch ip^' L ^ R carrying unit current and with group ve- 
locity directed towards region C (from C to oo). This is 

exemplified in Fig. 3 using the in mode (f>^J in . 

When the dispersive scale and the horizon surface grav- 
ity k are well separated (ui ~ n <C A) , the left-going mode 
does not significantly mix with the other three modes, all 
defined on the right-going branch of Eq. (5) [16], as we 
checked by numerical techniques based on [8]. Thus, the 
in- out scattering matrix is effectively a 3 x 3 matrix 



/ 



O- 1 



)-i 

— OJ 

( 2 ),i 

— OJ 



,(2) 



= A 



el 










( 2 ) 



t 



LU, OUt 



0-< 



(2), out 

— OJ 



(6) 

Given that the two ((p_ u )* have negative norms, the ma- 
trix coefficients satisfy normalizations conditions such as 

\<\ 2 -\^\ 2 -\^\ 2 = ^ (7) 

When working in the m-vacuum, the state without in- 
coming particles, the mean occupation numbers of outgo- 
ing particles with negative frequency are = |/3i^| 2 , 

whereas that with positive frequency is n" = + nf 2 ^, 
by energy conservation. That is, pair production here 
occurs through a two-channel Hawking-like mechanism. 
To approximately compute the coefficients of Eq. (6) in 
the regime k -C A we use connection formula techniques. 
We first decompose (j) u in both asymptotic regions L and 
R as a sum of plane waves: 

k = li + LSH^fy + L^{^ L y, 

The coefficients are connected by 



(8) 



(i) 

OJ 



TU 

-(1) 

( 2 ) 



(9) 
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where Ubh and L%h, respectively, describe the off-shell 
scattering on the two horizons [16] and [/hj describes 
the WKB propagation from one horizon to the other, 
i.e. it is diagonal and contains the exponential of iS® = 
i J dxk^(x), where is fc^, fc£, or k^. By off-shell we 
mean that these three matrices are not restricted to the 
two modes that govern the asymptotic scattering on each 
horizon considered separately, that is the growing mode 
is here kept in the mode mixing. In fact, since k^ has neg- 
ative imaginary part, e lS ^ is exponentially large. Simple 
WKB algebra shows that it grows as e AA where A is the 
distance between the two horizons. Concomitantly, since 
fcl — fct* ; e *s„ j g exponentially small. 

Let us show with an example how to determine the 



IS 



coefficients of Eq. (6). The globally defined 
constructed by imposing that the asymptotic amplitudes 
of the two incoming branches , (p_lj both vanish, 
see Fig. 3. Therefore, the three outgoing amplitudes 
are given by the second row of the matrix of Eq. (6). 
Moreover, using Eq. (9), these coefficients correspond to 

{RI,b£\b2 ] ) = {$\a£l,0) and (LJJ, , L^) = 
(0, 1, A_ u ). Solving the resulting system, we obtain 



/3 (1) = 
= -B 



bh x e ,5; 



BH 



x e 



' x a 



WH 



+ 0(e 4S i), 



x A 



WH 



0(e 



(10) 



where the a's and /3's in the above are the standard Bo- 
goliubov coefficients for black and white holes that en- 
code the thermal Hawking radiation [17]. By a similar 
analysis of other modes, all coefficients of Eq. (6) can 
be computed. Although the nonpositive definite conser- 
vation law of Eq. (7) does not bound these coefficients, 
the exponentially large factor in e AA cancels out from all 
of them. As a consequence, as can be seen in Eq. (10), 
the leading term is, up to some phase coming from the 
WKB propagation in region C, given by the Bogoliubov 
coefficients of Ubh and C/wh [16]. In the first two lines, 
one finds a product of two coefficients because the asso- 
ciated semiclassical trajectory passes through both hori- 
zons. Instead, in the third line only one coefficient is 
found because there is only a reflection on the black hori- 
zon. 

To identify the physical consequences of the pair cre- 
ation encoded in Eq. (6), we compute the expectation 
value of the stress-energy tensor 

_ 2 SS+ _ (0) (A) 

where T$) is the standard relativistic expression and 
arises from the Lorentz violating term of Eq. (3): 



~ (^ Q/3 <?W) 2 9^ 



(12) 



In the asymptotic region on the right of the white hori- 
zon, the field can be expanded as the superposition of 



the two right-going modes 



and 



,(l),out 



see Fig. 3, 



duj 4>[ 



l u.out * tA,OUt 



,(l),out~(l),out 



h.C. (13) 



In this region, the geometry is stationary and homoge- 
neous. Hence, the renormalized tensor T?j™ is obtained 
by normal ordering the above creation and destruction 
out operators. [The fact that there exist negative fre- 
quency asymptotic particles causes no problem in this 
respect. In fact, all asymptotic excitations have a pos- 
itive comoving frequency fl of Eq. (5).] Imposing that 
the initial state is vacuum, using Eq. (6), it is straight- 
forward to compute (0i n |T^™|0i n ) ■ The final expression 
contains an integral over lj of a sum of terms, each being 
the product of two modes 0£/ out , 0^'° ut and two coeffi- 
cients of Eq. (6). We do not need the exact expression 
because we only consider possible divergences. When 
lj > w max there are no negative frequency modes. Hence, 
the B coefficients of Eq. (6) vanish for lj > LJ max , and the 
stress-energy tensor cannot have ultraviolet divergences. 
Therefore, the only possible divergence can be found for 
lj -> 0. 

In each term of Tfi), there are two derivatives with 
respect to t or x, leading to two powers of lj, , or 



Analogously in T^t\ there are 4 powers of these. 



Now, from Fig. 2 we see that the wavenumbers k^\ k^ 
do not vanish for lj — > in L. Rather, they go to con- 
stant opposite values, that we call kg and — fc , respec- 
tively. Interestingly, these modes played a key role in [18] 
when studying the fluxes emitted by a white hole flow in 
an atomic Bose condensate. They are also associated 
with the macroscopic undulation observed in the experi- 
ments [19, 20] and theoretically described in Sec. III.B.3 
of [16]. Finally, the terms containing only spatial deriva- 
tives will not be suppressed for lj —> 0. The leading terms 
in (0i n |r^' ron |0i n ) are thus proportional to 



k 2 



4Trtt(k ) v g0 



dtu 



= (1) 



(14) 



The above integral gives the integrated mean occupation 
number of the two out species, and v g o is their asymp- 
totic group velocity in the t, X frame. The leading terms 
of (0i n |r/ 1 A ' ) ' rcn |0i n ) are proportional to Eq. (14) up to an 
extra factor of k 2 /A 2 . Since ko — A-^/uq — 1, one finds 
that the (0) and (A) components of the stress-energy ten- 
sor yield typically the same contribution. 



The key result comes from the fact that 



diverges 



as 1/lo 2 for lj — > 0, being the product of |/3 BH | 2 ~ 1/lj 
and |/3^ H | 2 ~ (This infrared behavior has been 

validated by numerical analysis.) Then, if the warp drive 
is created at some time, after some transient period the 
geometry will become stationary and the outgoing fluxes 
will reach their stationary values. However, only frequen- 
cies lj > 1/T will contribute after a lapse of time T. This 
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provides an infrared cutoff to the integral of Eq. (14), 
which gives an energy density scaling as 



/ duj 




ll/T 





That is, there is an infrared divergence that leads to a 
linear growth of the energy density. This result can be 
understood from the findings of [18]: the BH radiation 
emitted towards the WH horizon stimulates the latter as 
if a thermal distribution were initially present. In that 
case, it was also found that the observable (the density 
correlation function) increased linearly in T . 

Using quantum inequalities [21], it was argued [7] that 
k < 10 _2 ip , where tp is the Planck time, which implies 
that the growth rate is of the same order (unless A is 
very different from tp 1 ). In the presence of superluminal 
dispersion, warp drives are thus unstable on a short scale. 



III. SUBLUMINAL DISPERSION RELATION 

In this case, the warp drive system is also unstable. 
Indeed, the positive norm modes cftffl and (fiffl are now 
trapped between the two horizons in region C. Hence, 
they bounce back and forth and this induces an ex- 
ponentially growing amplification of the emitted radia- 
tion. As a result, the asymptotic fluxes, and therefore 
(0m|T/i™|0m), grow exponentially in time. This phe- 
nomenon is the subluminal version of the black-hole laser 
effect [22]. This dynamical instability is described by 
complex frequency eigenmodes that are asymptotically 
bounded [23]. In the original version, the analysis was 
performed with a superluminal dispersion relation and 
with a metric like that of Eq. (2) but with \V\ < 1 at 
±oo and |V^| > 1 in the central region (see added note). 
However, there is a precise symmetry between that case 
and the present one [16]. It basically consists in chang- 
ing both the sign of the dispersion relation and that of 
V + 1, i.e. interchanging sub and superluminal regions. 
This symmetry allows to infer that the complex frequen- 
cies governing the laser instability will share the same 
characteristics as in [23]. In particular, the growth rate 



of the most unstable mode, which is a non trivial function 
of the surface gravity of each horizon and their separation 
A, can be read of from the expressions of that work. 

IV. FINAL REMARKS 

In this paper, we generalized [7] by showing that super- 
luminal warp drives are unstable even when local Lorentz 
invariance is broken at very high energies. When the dis- 
persion relation is subluminal, the horizons act as a reso- 
nant cavity producing a dynamical black-hole laser insta- 
bility [22]. When it is superluminal, instead, the emitted 
flux grows linearly in time due to infrared effects. 

To conclude, we note that in [2] it was shown that 
close timelike curves can be obtained by combining sev- 
eral warp drives. Our results, together with those of [7], 
weaken that possibility because isolated warp drives are 
unstable irrespectively of the features of the dispersion 
relation in the ultraviolet regime. As a consequence, 
whereas former attempts to tackle the issue of chronology 
protection deeply relied on local Lorentz invariance [24], 
the present result suggests that this conjecture may be 
valid also for quantum field theories violating Lorentz 
invariance in the ultraviolet sector. 



Added note 

In fact, in the presence of superluminal dispersion, as 
it is the case in Hofava gravity [13], all black holes could 
be dynamically unstable. This conjecture is based on the 
observation (see Fig. 10 in [25] and the associated discus- 
sion) that any significant scattering in the region where 
\V\ > 1 induces signs of instability. It is therefore impor- 
tant to test this conjecture by a linear stability analysis 
based e.g. on the "improved" action of [26, 27]. 
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